Exponential increase of energy level density in atoms: Th and Th II 
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We present analytical estimates and numerical calculations showing that the energy level density 
in open-shell atoms increases exponentially with increase of excitation energy. As an example, we 
use the relativistic Hartree-Fock and configuration interaction methods to calculate the density of 
states of Th and Th II. The result is used to estimate the effect of electrons on the nuclear transition 
which is considered for the use in a nuclear clock. 
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I. INTRODUCTION 

The exponential increase of density of states is well 
known in nuclear physics (see, e.g. Similar problem 
in atomic physics has not been considered so far. It is 
well-known that the density of the Rydberg levels (with 
one excited electron) tends to infinity at ionization limit. 
However, there are also compound states with several 
excited electrons. As we will see below, the density of 
such states increases exponentially with the excitation 
energy. 

The concept of the density of states p(E) can be espe- 
cially useful for the cases of the dense and complicated 
spectrum. This takes place for practically any many- 
electron atom at sufficiently high excitation energies, in 
particular at energies close or above the ionization limit. 
In the latter case we speak about density of Feshbach 
resonances corresponding to the quasi-stationary com- 
pound states with several excited electrons. For atoms 
with open d or / shells the spectrum is complicated even 
at low energies. Experimental data for highly excited 
states are often absent. Accurate calculations using spe- 
cific atomic states in the area of dense spectrum is very 
difficult if not impossible. Even identification of such 
states is often problematic. On the other hand, statisti- 
cal approaches, which use functions obtained by averag- 
ing over large number of states, can be very useful. An 
interesting case was considered in Ref. Q where proper- 
ties of excited states of cerium atom were studied using 
statistical analysis. It has been demonstrated that the Ce 
atom is an example of quantum chaotic system in which 
highly excited states are similar to compound states of 
heavy nuclei. Density of states p is used in the criteria 
for the chaos to take place: Hij > Dij = 1/pij. Here 
is the off-diagonal matrix element of the effective Hamil- 
tonian between many-electron states i and j, and 
is the energy interval between the states which can be 
mixed by Hij . Statistical theory based on the properties 
of chaotic eigenstates @, Q allows one to calculate av- 
erage orbital occupation numbers, matrix elements and 
amplitudes between chaotic states, enhancement of weak 
interactions and increase of entropy. Another important 
example is the effect of atomic electrons on nuclear tran- 
sitions, which will be considered below. The concept of 
the density of states is also needed to estimate a number 



of states missed in experimental spectra. 

In this paper we present simple analytical estimates ex- 
plaining exponential increase of density of states contain- 
ing several excited electrons. Then we perform numerical 
calculations of density of compound states using thorium 
as an example due to its importance for the construction 
of nuclear clock. 

The proposed clock 0] utilizes the use of the nuclear 
transition between a metastable isomeric state and the 
ground state of 229 Th. The frequency of this transition 
is unusually low, about 7 eV, and can be accessed by 
lasers. Additional interest is due to the fact that this 
transition must be very sensitive to the time-variation 
of the fine structure constant [f|. The effect of elec- 
trons on the nuclear transition in Th II and Th IV ions, 
also known as electronic bridge process, was considered 
in Refs. [1,0] (see also [§[). In this process the nuclear 
transition is accompanied by an electron transition from 
the ground state to an exited state followed by emis- 
sion of the photon. The electronic bridge process may 
have significant effect on the performance of the nuclear 
clock in the case of resonance between nuclear and atomic 
transitions. However, neither nuclear frequency nor elec- 
tron spectrum are known sufficiently well to make reli- 
able estimations. The measured values for the frequency 
of the nuclear transition varies form 3.5 ±1.0 eV @ to 
7.6 ± 0.5 [l3|. Electron spectrum of Th II is also not 
known at energies around 7.6 eV. In this situation instead 
of calculating the value of the electronic bridge process 
one can consider estimation of the probability for the re- 
sult to be within some range of values. This can be done 
if density of electron states as a function of the energy 
is known. Knowledge of the electron spectrum density 
in Th II is also needed for the accurate measurements of 
the nuclear transition frequency. E. Peik and collabora- 
tors [ll[ plan to use laser excitation of a Th II electron 
state which may transfer part of the excitation energy to 
induce the nuclear transition. 



II. SIMPLE ANALYTICAL ESTIMATE OF THE 
NUMBER OF COMPOUND STATES 

Let us consider an open ground electron shell con- 
taining n electrons and go single-electron states; usually 
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go = J^i 2(2/ + 1) where I is the orbital angular momen- 
tum of a single-electron state within this shell. A total 
number of many-body states within this shell is equal to 
the binomial coefficient, 



90 n\(go — n)\ 



cxp [nln (ffp/n) + 1] 
V(27rn) 



(1) 



In the last expression we assumed go S> n and used the 
Stirling formula 



(2) 



which gives an excellent accuracy even for n = 1 (the 
correction is — l/(12n)). Assume for simplicity that all 
these many-body states have approximately the same en- 
ergy. Let us now consider next electron shell which is 
separated from the ground state shell by a single-particle 
energy interval ui and has g\ single-electron states . If a 
total excitation energy of a many-electron state is equal 
to E, we may transfer k = E/w electrons from the ground 
shell to the next shell. Now we have a total number of 
many-electron states 



N(E) = C%> + CJ!C a n °_ k C^ . 



(3) 



Here the first factor C% gives the number of ways to 
select k excited electrons from the available n electrons, 
the second factor gives the number of many-body states 
in the ground shell and the last factor gives the number 
of many-body states in the next shell. For n^> k 



N(E) 



N \1 + 

No + | exp (bE) . 



exp [k In (n 2 gi/k 2 go) + 2k] 
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where 



cxp [nln (g /n) + 1] 



ln[(n 2 ffl )/(fc 2 5o )] 



(4) 



(5) 



(6) 



This rough estimate shows that the number of compound 
states increases exponentially with the number of elec- 
trons in the open shell n and with the number of excited 
electrons k = E/uj. Inclusion of the electron excitations 
to higher shells and the Rydberg spectrum of the single- 
electron excitations will only make the number of states 
larger. 

The derivation of formula ((4]) assumes that the number 
of electrons n in the open shell is large. We will see below 
that the exponential increase of the density of states takes 
place even for n as small as three. 



TABLE I: Number of states of even parity and total momen- 
tum J=1.5 in 10000 cm" 1 energy intervals in Th II. 
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Exp." 


Full CI 6 
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±1aa 
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10000 


5 


5 


1 


3 


10001 


20000 


3 


3 


5 


5 


20001 


30000 


4 


5 


4 


7 


30001 


40000 


10 


11 


7 


15 


40001 


50000 


18 


10 


19 


8 


50001 


60000 


8 


18 


5 


26 


60001 


70000 




50 


44 


90 


70001 


80000 




93 


94 


160 


80001 


90000 




161 


173 


300 


90001 


100000 




128 


142 


168 



"Experiment [1311 
6 Short basis 

c No diagonalization, short basis 
d No diagonalization, long basis 



III. NUMERICAL CALCULATIONS 

Density of atomic states p can be defined via the rela- 
tion 



N(E) = / p(e)de, 
Jo 



(7) 



where N(E) is the number of atomic states in the interval 
from zero energy, which corresponds to the ground state, 
to the energy E. Total density p is the sum of the partial 
densities p p which correspond to the states of definite 
total momentum J and parity: 



(8) 



We use relativistic Hartree-Fock (RHF) and configura- 
tion interaction (CI) methods to calculate the number 
of states for given J and parity. The partial density is 
calculated by the numerical differentiation 



Pp(e) 



N p (e + S) -N p (e-S) 
25 



(9) 



The value of 6 is chosen to have a smooth function for 
p p . This is achieved when the number of states in the 
numerator is large. The density function defined this 
way represents inversed energy interval D p averaged over 
large number of neighboring states, p p — 1/D p . 

We perform the configuration interaction (CI) calcu- 
lations using the V^ -1 basis (see e.g. [Hj]) calculated 
in the frozen RHF field of the [Ra]6c? 2 configuration for 
Th II and the [Ra]6d 2 7s configuration for Th I. To test 
the accuracy of the calculation we use two different basis 
sets for Th II: a short one and a long one. The purpose of 
using short basis is to show that the non-diagonal matrix 
elements of the CI matrix can be neglected in the evalu- 
ation of the density of states. This allows us to use much 
larger basis and, consequently, move to higher energies. 



3 



In the short basis we include six lowest RHF states 
above the core in each of the partial waves up to l max = 4. 
Two more states in each partial wave are added in the 
long basis. Apart from that, we allow only single and 
double excitations from a reference configuration to con- 
struct three-electron states for the CI calculations with 
the use of the short basis. Three-electron excitations are 
allowed when long basis is used. 

For the short basis we perform both the full-scale CI 
calculations which include matrix diagonalization, and 
the simplified calculations in which off-diagonal matrix 
elements of the CI matrix are neglected. This corre- 
sponds to an approximation Em ~ (M\H CI \M), where 
\M) is a three-electron CI state of a definite total momen- 
tum J and parity, H CI is the CI Hamiltonian. For the 
long basis we neglect the off-diagonal matrix elements. 
The use of this approximation allows us to deal with large 
number of states (~ 10 5 ). On the other hand, neglecting 
the non-diagonal CI matrix elements has little effect on 
the density of states. 

Table U shows experimental and calculated number of 
states in the 10000 cm" 1 energy intervals of Th II. We 
consider even states of total momentum J=1.5 as an ex- 
ample. These states are important for the electronic 
bridge process |7j. Experimental data ends at about 
56000 cm" 1 . Next two columns correspond to exactly 
the same CI matrix obtained with the use of the short 
basis set. It has been diagonalised in the full CI calcula- 
tions while only diagonal matrix elements are considered 
in the next column. One can see that neglecting the off- 
diagonal CI matrix elements does not lead to dramatic 
changes in the density of states at high energies where the 
density function is well defined (see eq. © and discussion 
below it). 

The last column of Table Q] shows density of states ob- 
tained with the use of long basis without diagonalization. 
The use of long basis adds more states at large energy 
and has little effect on the number of states at low energy. 
Therefore, we use only long basis in the calculations for 
Th I. 



IV. RESULTS 

The results of the calculations for the total density of 
states are shown on Fig. [TJ It is also shown that the 
calculated data is fitted very well by simple exponential 
functions. We use two fitting formulae 

p{e) = - exp(e/e ) (10) 

and 

p{e) = — exp(e/e ), (11) 
e 

where A and eo are fitting parameters. Formula (|11[) 
is more consistent with the analytical formula Q while 




Excitation energy (a.u.) 

FIG. 1: Total density of states of Th I and Th II as func- 
tions of excitation energy. Points show the result of the CI 
calculations, solid line - eq. QD) , dashed line - eq. The 
density in Th I is higher since the number of electrons in 
the open shell is larger, n = 4, and the interval between the 
single-particle levels ui is smaller - see eq. Q. 



TABLE II: Fitting parameters A and eo (a.u.) for total den- 
sities of states of Th II and Th I. 





Formula (Tl0|) 




Formula (TiT]) 




A e 


A 




Th II 


9.8 0.0625 


6 


0.045 


Th I 


1.06 0.025 


2 


0.021 



formula (JTTJJ) is simpler and behaves better at e 0. 
The values of the fitting parameters A and eo for the 
total density of states of Th II and Th I are presented in 
Table ITI1 It turns out that eo can be kept the same for all 
partial densities. The values of parameter A p for partial 
densities are presented in Table Mil 

The density of states must go to infinity when energy 
is approaching ionization potential. This is due to the 
infinite number of Rydberg states which all have about 
the same energy, close to the ionization potential. In 
contrast, the number of calculated states goes down at 
large energies. This is the cut-off effect of the single- 
electron basis. Adding more highly excited orbitals into 
the basis increases the number of many-electron states at 
high energy (see Table [TJ) . Ionization potential of Th I is 
50867 cm" 1 = 0.23 a.u. Q3 while the one for Th II is 
about 96000 cm" 1 = 0.44 a.u. [l5j]. As can be seen from 
Fig. [TJ the fitting works very well for the energies very 
close to the ionization potentials of Th II and Th I. At 
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TABLE III: Fitting parameters A (formula (|10[0 for partial densities of states of Th II and Th I. 











Even states 












Odd states 






Total 


Th II 


J 


0.5 


1.5 


2.5 3.5 


4.5 


5.5 




0.5 


1.5 


2.5 


3.5 


4.5 


5.5 








A 


0.6 


0.9 


1.2 1.2 


0.9 


0.7 




0.4 


0.8 


0.9 


0.9 


0.8 


0.5 




9.8 


Th I 


J 





1 


2 3 


4 


5 


6 





1 


2 


3 


4 


5 


6 


7 




A 


0.03 


0.07 


0.08 0.07 


0.07 


0.05 


0.02 


0.03 


0.09 


0.13 


0.14 


0.12 


0.09 


0.05 


0.02 1.06 



higher energies more basis states must be included. 

Now we can use the density function to estimate the 
probability of the electronic bridge process in Th II. Ac- 
cording to Ref. Q the ratio /3 of the electron bridge width 
to the width of the direct nuclear transition (the electron 
enhancement factor) is 

B*( —Y ^ (12) 

\uj n J 3(2J n + l)(2Ji + l){u in + uj N f' V ' 

where n is the resonance electron state, uj is the photon 
frequency , ujn is the frequency of the nuclear transi- 
tion, J„ is the total angular momentum of the resonance 
state, Ji is the total angular momentum of the ground 
state (J n = Ji=1.5 in our case), R n is the combination 
of the electron matrix elements for electric dipole and 
hyperfine transitions between atomic ground state, reso- 
nance state n and some odd-parity low lying final atomic 
state, uji n is the frequency of the electron transition from 
the ground to the resonance state, ujn is the frequency 
of the nuclear transition. Resonance corresponds to the 
condition uji n — wjv- As it was discussed in the intro- 
duction the exact value for the nuclear frequency is not 
known and the experimental values range from 3 eV to 
8 eV. Estimation of the effect for 3.5 and 5.5 eV were 



performed in Ref. |7j. The effect for higher con was not 
presented since the electron spectra are not known for 
such energies. Here we perform the estimation for the 
latest measured value ujn — 7.6 eV = 0.28 a.u. Using 
formula (fTU|) and the data from Table Mil we find that 
the average energy interval between the relevant electron 
states is D p — l/p p — 6 x 1CT 4 a.u at this energy. The 
energy denominator in (TT2l is at least two times smaller. 
Therefore, the electron enhancement factor is 

/5~2x 10 5 i?„ « 4 x 10 3 . (13) 

Here we use the root mean square value for R n which is 
R n s» 2 x 1CU 2 a.u. for the energies above 7 eV 0j- 

An interesting effect may take place for neutral tho- 
rium if the frequency of the nuclear transition is smaller 
than the ionization potential of the atom. The density 
of states of Th I is huge and the resonance situation is 
highly probable. This can be used to excite the nucleus in 
a two-step process in which a resonance electron state is 
excited first, then it gives its energy to excite the nucleus. 
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